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Cluster abundance measurements are among the most sensitive probes of the amplitude of matter 
fluctuations in the universe, which in turn can help constrain other cosmological parameters, like 
the dark energy equation of state or neutrino mass. However, difficulties in calibrating the relation 
between the cluster observable and halo mass, and the lack of completeness information, make this 
technique particularly susceptible to systematic errors. Here we argue that a cluster abundance 
analysis using statistical weak lensing on the stacked clusters leads to a robust lower limit on 
the amplitude of fluctuations. The method compares the average weak lensing signal measured 
around the whole cluster sample to a theoretical prediction, assuming that the clusters occupy the 
centers of all of the most massive halos above some minimum mass threshold. If the amplitude 
of fluctuations is below a certain limiting value, there are too few massive clusters in this model 
and the theoretical prediction falls below the observations. Since any effects that modify the model 
assumptions can only further decrease the prediction of the model in the context of this method, the 
limiting amplitude becomes a robust lower limit. Here, we apply it to a volume limited sample of 
16,000 group/cluster candidates identifled from isolated luminous red galaxies (LRGs) in the Sloan 
Digital Sky Survey (SDSS). We flnd 0-8(^^/0. 25)°'^ > 0.62 at the 95% c.l. after taking into account 
observational errors in the lensing analysis. While this is a relatively weak constraint, both the 
scatter in the LRG luminosity-halo mass relation and the lensing errors are large. The constraints 
could improve considerably in the future with more sophisticated cluster identification algorithms 
and smaller errors in the lensing analysis. We argue that the existence of a lower limit from cluster 
abundance is rather general, and demonstrate that Malmquist bias dominates over Eddington bias 
in this type of analyses. 



PACS numbers: 

I. INTRODUCTION 

The abundance of clusters of galaxies has been long 
recognized as potentially one of the most powerful probes 
of cosmological parameters. The main idea is that the 
cluster abundance can be related to the abundance of 
dark matter halos, which are compared to the theoreti- 
cal halo mass function predictions. These have an expo- 
nential cutoff at the high mass end above the so-called 
nonlinear mass. The nonlinear mass is most sensitive 
to the amplitude of fluctuations, usually expressed as 
(Tg, but also depends on other cosmological parameters 
such as the matter density Qm- By measuring the clus- 
ter abundance evolution with redshift, one can determine 
the growth rate of structure. Combination with other 
probes such as the cosmic microwave background, super- 
novae, and galaxy clustering makes cluster abundance 
measurement a very powerful probe of such fundamen- 
tal parameters as neutrino mass and dark energy density 

Qiii. 

Existing cluster surveys identify clusters based on an 
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observable property like luminosity. One typically identi- 
fies a flux or surface brightness limited sample in a given 
area based on one of the observables, and follow up obser- 
vations determine the cluster redshifts via cluster mem- 
ber spectroscopy. With redshift information, volume- 
limited samples of constant luminosity threshold can be 
defined, and the cluster volume density determined. In 
X-ray surveys, the total cluster luminosity is used as an 
observable [5j, although other variables may reduce the 
scatter Q. In optical surveys, either cluster richness, to- 
tal luminosity or central galaxy luminosity is used as the 
cluster observable 0, H, @- Finally, Sunyaev-Zel'dovich 
(SZ, 10]) surveys use either the total SZ luminosity or 
the central decrement. To relate the cluster sample to 
the halo mass function, one must connect the observable 
to the halo mass. With X-ray data, the halo mass can be 
determined from the measurement of X-ray temperature 
and intensity profile, which allows one to determine the 
mass profile through the equation of hydrostatic equilib- 
rium. Optical surveys may use velocity dispersion mea- 
surements, but recently weak lensing calibration has also 
become possible [ll|. The SZ signal may already be a 
reliable tracer of the dark matter halo mass, since it is 
only weakly affected by astrophysical complications like 
gas cooling, feedback, cosmic rays, etc. [H, [3]. While 
this technique is still in its infancy, several surveys should 
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transform it into reality in the near future. 

In all cosniological applications of cluster abundances, 
an accurate calibration between the observable and dark 
matter mass is required. It is important not just to de- 
termine the mean of this relation (e.g. [13])) but also its 
scatter [l^ EE] and higher moments. Small errors in this 
calibration may significantly affect the final cosmological 
parameters. For example, given the steepness of the mass 
function, a Gaussian scatter brings many more low mass 
halos into the cluster sample than large mass halos out 
of the sample, increasing the abundance for a fixed cutoff 
in the observable. On the other hand, if a cluster is truly 
dark, it cannot be observed at all, so the cluster sample is 
incomplete at a given mass threshold, and the observed 
abundance is below the true value. In general, given the 
steepness of the mass function, even a small deviation 
from Gaussianity at the tail of the error distribution can 
cause a significant effect, which is very difficult to identify 
with current methods using small subsamples of clusters 
for which the mass is determined individually. 

The second problem in cluster abundance analyses is 
that some mass determination methods may be unreli- 
able because the assumptions on which they rely may be 
violated. For example, mass determination from X-ray 
measurements of gas intensity and temperature using hy- 
drostatic equilibrium assumes that the pressure support 
is thermal, while additional sources of pressure (cosmic 
rays, magnetic fields, bulk motions, etc.) are ignored 
|13.[l7|. Other X-ray observables (e.g. the product of the 
X-ray temperature and the gas mass within rsoo, or Y^, 
) may not require the assumption of hydrostatic equi- 
librium, but still must be calibrated against simulations, 
which may not include all of the relevant physics. It is 
generally accepted that gravitational lensing is the most 
direct tracer of the halo mass, so any attempts at deriv- 
ing reliable constraints from cluster abundances should 
ideally be based on lensing. 

The third problem is that the calibration between the 
observable and the halo mass is usually done on a sub- 
sample of clusters, which may not be representative. An 
example is the relation between mass and X-ray lumi- 
nosity, which is often calibrated on relaxed clusters for 
which the hydrostatic equlibrium assumption may be 
valid, while the results are then applied to the whole 
cluster sample even though non-relaxed clusters may not 
obey this relation. 

These arguments suggest that only if clusters are iden- 
tified by mass selection, such as when using the dark 
matter maps derived from gravitational lensing, can one 
avoid these potential systematic errors. Unfortunately, 
mass selection of clusters is still far away from being 
practically implemented. Also, while in principle a mass- 
selected sample lacks the difficulties of the current meth- 
ods, the amount of cosmological information is limited 
without redshift information. If the redshift is sought af- 
terwards, there is again the possibility that the selected 
cluster candidate lacks bright galaxies, bringing in astro- 
physical complications. Moreover, there is considerable 



scatter between the halo mass derived from lensing and 
the more traditional halo mass estimates jl^. For ex- 
ample, in some cases this mass selection may result in a 
random superpositions of filaments rather than true viri- 
alized clusters, so one must be able to quantify the rate 
of such occurrences. Alternatively, several promising self- 
calibration methods were proposed for future large data 
sets [l^ [131 , and in some cases were already applied to 
existing data sets fi\ , but it remains to be seen whether 
these methods can reliably solve all of the aforementioned 
problems. 

It is worth asking if there is any robust information 
that can be extracted from the cluster abundances with- 
out making additional assumptions that increase the sus- 
ceptability to hidden systematic errors. We have estab- 
lished above that to derive a robust constraint, we must 
use gravitational lensing to relate the observable to mass; 
that this analysis should use the complete sample of clus- 
ters; and that the effects of scatter and incompleteness 
are difficult to fully account for. Here, we propose a 
method that uses lensing analysis on the entire cluster 
sample. Because the issues of completeness and scatter 
cannot be fully solved without additional assumptions, 
only an inequality can be established. 

We propose the following analysis. Given a cluster 
sample with a known number density at a given redshift, 
we compute the halo mass function for a given cosmol- 
ogy (i.e., amplitude of fluctuations assuming a value for 
ri„i and other cosmological parameters). From this, we 
can compute the theoretical prediction for the mean weak 
lensing signal around the clusters assuming that the clus- 
ter positions are at the centers of all of the most massive 
halos predicted by the model (with a minimum halo mass 
Mfnin determined by the cluster abundance) . This is the 
maximally allowed lensing signal for that cosmological 
model, since placing the clusters at any other position 
can only lower the lensing prediction. From the data, 
one can compute the mean weak lensing signal around 
the clusters by stacking the individual signals. Stack- 
ing reduces the noise, so that the resulting S/N can be 
high even if it is low from an individual cluster. Further- 
more, it avoids the problems with individual cluster lens- 
ing mass determination mentioned in fl§|, since we only 
require that the stacked signal give information about 
the cluster mass in the mean. As in (2l| . comparison 
with N-body simulations suggests that this assumption 
is valid. 

This signal can then be compared to the theoretical 
predictions. These depend mostly on the amplitude of 
fluctuations: lowering it will reduce the number of high 
mass clusters and the lensing signal amplitude. At some 
point the theoretical predictions fall below the observa- 
tions. Such a cosmological model cannot be allowed, 
since any failure of the assumption that the cluster sam- 
ple corresponds to the centers of the most massive ha- 
los can only reduce the lensing signal, making the dis- 
crepancy worse. The amplitude of fluctuations for which 
the theoretical prediction still matches the observations 
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within the observational errors is thus a robust lower 
limit. 

This method can be applied to any cluster sample, no 
matter how incomplete, but the derived lower limit im- 
proves if the sample is closer to the complete sample of 
most massive halos. This may be a small effect if one 
is on the exponential tail, where a small change in am- 
plitude causes a large change in the abundance. Note 
that within the context of this method, we cannot estab- 
lish when the inequality becomes equality. This would 
require an assessment of completeness and scatter in the 
observable-mass relation, which as argued above is very 
difficult and may not even be possible if some of the clus- 
ters are truly dark. 

In this paper, we apply the method to a cluster sam- 
ple derived from isolated luminous red galaxies (LRGs) 
in Sloan Digital Sky Survey (SDSS), for which the high 
signal to noise of the lensing signal has already been es- 
tablished . The advantages of this sample are that it 
covers an unprecedently large volume and is essentially 
volume limited. On the other hand, the LRG luminosity 
is unlikely to be a perfect tracer of cluster mass, so with 
a better tracer the constraints may be improved in the 
future (for example using the maxBCG cluster sample 
presented in (Hi)- 

Here we note the cosmological model and units used 
in this paper. All data-related computations assume 
a flat ACDM universe, though Sl^ itself is allowed to 
vary. Distances quoted for transverse lens-source separa- 
tion are comoving (rather than physical) /i~^kpc, where 
Hq = 100 /i kms~^ Mpc^^. Likewise, AE is computed 
using the expression for Ej^ in comoving coordinates, 
where Ec is the critical surface density that depends on 
the distance ratios [0]. For the typical lens and source 
redshifts in this study, Ec depends on cosmology only at 
the 0.1% level even for extreme changes in cosmology. 
In the distance units used, Hq scales out of everything, 
so our results are independent of this quantity. When 
masses are quoted in Mq rather than H^^Mq, we have 
used h = 0.7. 



II. DATA SAMPLE AND ANALYSIS 

As the basis of the cluster sample we use the spectro- 
scopic sample of LRGs from the latest release of SDSS, 
data release 5 [H, |23| ■ These are good tracers of mas- 
sive halos, as suggested by their large bias [l^. Previous 
lensing analyses [9| have established that the minimum 
halo mass is ~ 5 x IO^^Mq and that there is a correla- 
tion between LRG luminosity and halo mass, M (x L^, 
for masses up to 2 x 10^** M0. The complete, approxi- 
mately volume-limited sample consists of ^40,000 galax- 
ies with 0.15 < z < 0.35 and a comoving volume of 
0.44(Gpc//i)'^. We identify a higher-luminosity subsam- 
ple of host LRGs using spectroscopic galaxy counts in 
cylinders of comoving radius 1 /i~^Mpc and line-of-sight 
length Av — ±1200 km s~^. Since we only want to elimi- 



nate galaxies for which there is another, brighter spectro- 
scopic LRG nearby, we simply require that the LRGs in 
our sample either (a) be the only one in the cylinder cen- 
tered on its position, or (b) be the brightest in the cylin- 
der. This cut, which is weaker than the overly-stringent 
cut in using a 2 h~^Mpc cylinder radius, eliminates 
about 7.5% (rather than 13%) of the complete sample. 

For the LRG luminosity, we use r-band model mag- 
nitudes k + e-corrected to 2; = as in to cre- 
ate luminosity-threshold samples with Mr < —22.3, 
< -22.6, and < -23 containing 15 635, 5 099, and 902 
galaxies (respectively). Because of scatter in M(L), there 
is no guarantee that this will have the desirable effect of 
selecting the most massive among the halos in the sam- 
ple. 

LRGs targeted for deeper spectroscopy in SDSS have 
a very specific color selection [26| , and it is unlikely that 
only the galaxies that satisfy the LRG criteria can be 
central galaxies in clusters above our mass threshold. In- 
deed, in m, only 70% of SDSS maxBCG clusters in re- 
gions with spectroscopy contain LRGs with spectra. We 
quantify the effects of color-related incompleteness us- 
ing a comparison against the SDSS Main galaxy sample, 
which is purely flux-limited at a brighter limiting flux. 
The highest luminosity sample from our previous Main 
sample lensing analyses (27j , with Petrosian r-band mag- 
nitude /c-corrected to z = 0.1 satisfying —22 > Mr > —23 
(the L6 sample), overlaps significantly with the spectro- 
scopic LRG sample in redshift, and thus is an ideal sam- 
ple for quantifying the effects of the LRG color cut. LRGs 
constitute 65% of this brightest Main subsample, so one 
possibility is to supplement the LRGs with the galax- 
ies that are equally as bright but not as red. The Main 
sample limiting apparent magnitude is 17.77, so these L6 
galaxies appear significantly brighter than the complete 
LRG sample (r < 19.1), partially because their maximum 
redshift is lower (~ 0.25 rather than 0.35). However, 
their lensing-weighted mean redshift, (z) ~ 0.2, is not 
that different from that of the LRG sample, (z) ~ 0.24, 
because lensing gives the highest signal for the galaxies 
that are halfway between z = and the sources at red- 
shift 0.3-0.6. The L6 Main sample galaxies were split into 
the 65% that pass the LRG selection and the remaining 
35% that fails it. 

As shown in figure [U weak lensing analysis reveals 
that halos of the non-LRG L6 galaxies are (30 ± 20)% 
less massive on average. The masses were derived from 
NFW ([11]) fits to the inner 0.5/i-iMpc with fixed con- 
centration c = 6, where host group/cluster contributions 
for satellite lenses are negligible (since we do not have 
a reliable environment estimator for L6 that will allow 
us to remove satellite galaxies as we do for the LRGs). 
The reduced values for these fits are 0.53 (non-LRGs) 
and 1.0 (LRGs), which suggests that the NFW model is 
an acceptable description of the data (particularly when 
accounting for the increase in due to noise in the jack- 
knife covariance matrices, (3g|). Unlike in Q, the virial 
radius in this case is defined as the radius within which 
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FIG. 1: The lensing signal for the brightest Main sample 
galaxies, divided into those that are spectroscopic LRGs and 
those that are not. We also show the best-fit lensing signal 
on small scales (where the contribution due to host halos for 
those that are satellites is small), giving average estimated 
virial halo masses using the methodology of Q. 



FIG. 2: The lensing signal for the spectroscopic LRGs in lu- 
minosity threshold samples, Mr < —22.3 (black hexagons), 
< —22.6 (red crosses), and < —23 (blue circles). Bootstrap 
la errors are shown, along with the theoretical signal for the 
best-fit (Tg with fim ~ 0.25 (fixed) as a solid line, and the 95% 
CL lower limit as a dotted line. 



the average halo density is 280p. Allowing the concen- 
tration to vary lessens the difference to ~ 10%, as the 
best-fit concentration for the L6 LRG sample is higher 
than that of the L6 non-LRG sample. It is important 
to note that these non-LRG L6 galaxies are on average 
~ 0.1 mag fainter than L6 LRGs. Given that we have 
previously established M cx at the high-luminosity 
end P, we would indeed expect the L6 non-LRGs to be 
20% less massive. Within the errors, these results there- 
fore justify our assumption that non-LRG galaxies are in 
equally massive halos as LRG galaxies of the same lu- 
minosity in r. Thus, we increase the number density of 
LRGs by 38% for Mr < -22.3, and 25% for < -22.6 
and Mr < —23 according to the fractions of L6 galaxies 
that are classified as LRGs in these bins. We use these 
increased abundances (including the correction for the 
overly-stringent satellite cut described previously) in this 
paper while analyzing the LRG sample signal presented 
in |9| . The abundance corrections due to the satellite cut 
vary based on LRG luminosity, and lead to an increase of 
7.5% in the abundance of the complete sample, or 6.0%, 
2.6%, and 0.7% for the three luminosity threshold sam- 
ples. 

The weak lensing analysis is the same as in |9|. More 
than 30 million source galaxies are identified, their shape 
measurements obtained using the Reglens pipeline, in- 
cluding PSF correction done via re-Gaussianization [2^ 
and with cuts designed to avoid various shear calibra- 
tion biases. A full description of this pipeline can be 



found in [30|. The main difference relative to that pa- 
per is that the results of the Shear TEsting Programme 
(STEP) comparison to simulations are now available and 
suggest that Reglens is well calibrated at the 2% level 
[31]. Our redshift distributions are calibrated using the 
DEEP2 spectroscopic survey [13, HI]- As in [30], we as- 
sume a total 8% calibration uncertainty (Icr), though this 
is likely to be a conservative estimate. 

Since we know the LRG lens redshifts, we can express 
the lensing signal in terms of the differential surface mass 
density AE as a function of transverse separation R. We 
average the signal over all the clusters, so the result is the 
average weak lensing profile around them. Our lensing 
signal estimator, with all associated tests and corrections 
(for imperfect sky subtraction, intrinsic alignments, non- 
weak shear, and other effects), is described in detail in 
The results are shown in figure [2] for various luminosity 
threshold subsamples. We see that the signal increases 
with luminosity, so the luminosity thresholds do have the 
desirable effect of preferentially selecting more massive 
halos. We also see that the signal is inconsistent with zero 
at several standard deviations for all luminosity threshold 
samples and transverse separations shown in the figure. 
The total S/N (over all transverse separations using the 
full jackknife covariance matrix) is 19, 17, and 11 for the 
Mr < -22.3, < -22.6, and < -23 threshold samples 
respectively. 
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III. RESULTS 

We want to compare the weak lensing signal to a theo- 
retical prediction under the assumption that all of the 
most massive halos have been occupied. We use the 
Sheth and Tormen halo mass function [s^, but with pa- 
rameters as suggested in fs^l and defining the halo mass 
using the radius within which the average density is 280p. 
These two choices represent departures from the analy- 
sis in 0. The halo mass function specifies the number 
density of halos as a function of mass, dn/dM. For each 
luminosity threshold we select the minimum mass such 
that the integrated number density equals the observed 
one, (dn/dM)dM — uq, where uq is 5.20, 1.48, and 

0.26 X 10~^(/i/Mpc)^ for the three luminosity threshold 
subsamples (with = 0.25). 

Once we have the minimum mass, we compute the pre- 
diction for the weak lensing signal, 

where AI](Af) is the differential surface density for a 
cluster of mass M, which we model as an NFW profile 
with concentration parameter that depends on the non- 
linear mass and which must be calculated for each cos- 
mology separately using c(M) = 10{M/Mni)~°-'^^ [Si]. 
Here Mni is the nonlinear mass, defined as the mass in 
a sphere within which the rms overdensity fluctuation is 
1.68. The concentration dependence on the cosmological 
model mildly enhances the scaling with amplitude, be- 
cause in a model with a lower amplitude of fluctuations, 
not only will Mmin at a given abundance be lower, but so 
will be concentration, which reduces the signal at scales 
smaller than the virial radius. As in our previous analy- 
sis with this sample, we include a stellar component with 
mass and radius that is fixed to a value determined from 
the optical luminosity. This component is important for 
0.04 < R < 0.08/i~^Mpc, a small fraction of the region 
used for the fits (0.04 <R< 2h-^Mpc). 

This calculation is done separately for each cosmology. 
For the mass function, the most important parameter is 
the nonlinear mass, which depends on the amplitude of 
fiuctuations as and mass density flm- For erg — 0.75, 
ilm — 0.25, z — 0.24, Us — 1, and h = 0.7, figure [3] shows 
the effect of changing these cosmological parameters on 
our chosen mass function and confirms the dominant im- 
portance of CTg and Qm- The shape of the power spectrum 
also enters the mass function calculation [s^, but the 
variation in the mass function (figure[3]) is small given the 
current uncertainties on the shape of the power spectrum. 
For example, the change in the mass function is three 
times smaller for Ug than for when using WMAP3, 
supernovae, galaxy clustering, and Lyman-a forest, ^t\ . 
Our use of ris = 1 gives conservative bounds on cg, in the 
sense that decreasing it by 4% to the value from WMAP3 
would decrease the expected abundances by that amount, 
increasing the best-fit CTg by ^ 2%. Because it is approxi- 
mately degenerate with ilm for M < 10^'^h~^MQ, for the 
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FIG. 3: Parameter dependence of the mass function at the 
LRG mean redshift for the fiducial cosmology described in 
the text. 



fainter threshold samples where those halo masses dom- 
inate, one can estimate the best-fit as for a 4% lower 
ng by simply lowering the fiducial fl„i by that amount 
instead. 

Since we fix /i = 0.7, we effectively have a different 
power spectrum shape for each value of f2m, but as shown 
in the figure, the effect of Qm on the nonlinear mass 
definition is the main source of our sensitivity to this 
parameter, so the curves of variable f2„i with and without 
fixing the shape parameter F — fim/i are similar. Since 
the matter density affects a number of quantities such as 
the calculation of distances and volumes, the change in 
growth factor from z = 0.24 to z = 0, and the definition 
of virial radius, we perform the analysis separately for 
three values, ilm = 0.20, 0.25, 0.30. For each value of ilm, 
we vary only as while fixing the remaining parameters. 

We perform the analysis using minimization with 
respect to as- If as is high, then we will have many high 
mass halos, so the predicted lensing signal will exceed the 
observed one. Since deviation from the assumption that 
all of the most massive halos are in our sample can only 
reduce the signal, such a model cannot be excluded. Low- 
ering the amplitude of fluctuations reduces the number 
density of massive halos, so the lensing signal prediction 
decreases. As a result, there is a value that fits the data 
best. Reducing the amplitude further drops the model 
predictions below the observations, so such a model can- 
not be resurrected. Figure [2] shows the results of this 
fitting procedure, showing that we can obtain a good fit 
for each luminosity threshold. The best-fit reduced 
values for each luminosity threshold are 1.4, 1.2, and 1.2 
from faintest to brightest. Accounting for the noise in the 
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covariance matrix, these values give a p{> x^) (prob- 
ability to exceed the value by chance if the model is a 
good description of the data) of 20%, 40%, and 40%. We 
also show the prediction of a model which is excluded at 
the 95% confidence level. 

To determine errors on the lensing signal, we divide the 
survey area into 200 bootstrap subregions, and generate 
2500 bootstrap-resampled datasets. We repeat the anal- 
ysis for each resampled dataset, introducing the 8% sys- 
tematic calibration uncertainty at this stage. The final 
outcome of our analysis is a probability distribution for 
(jg for each luminosity threshold. This procedure avoids 
the problems [ssl [soj of overly-optimistic parameter con- 
straints due to noise in the bootstrap covariance matrix 
which tends to increase the fit so it deviates from the 
usual distribution, and leads to broader parameter dis- 
tributions than we would have obtained by naively using 
Ax^ values to find confidence regions. The distributions 
are plotted on figure [3] for = 0.25. We see that the 
luminosity thresholds Mr < —22.3 and < —22.6 give the 
strongest constraints on ag, while for < —23 the con- 
straints weaken. If we had a complete sample with no 
scatter then we should get the same constraints for all 
the samples (assuming the effects from cosmological pa- 
rameters besides amplitude are negligible), so the weaker 
constraints at the bright end suggest that the relation be- 
tween LRG luminosity and halo mass breaks down there. 

We then repeat the analysis for other values of fim, 
covering the range between 0.2 and 0.3. We find that 
we can match the probability distributions by intro- 
ducing a variable (T8(ilm/0.25)°''^. The final outcome 
from our analysis is thus the probability distribution on 
this parameter. The constraints can be expressed as 
erg (^^m/0. 25)"-^ > 0.68 at 50% c.L, crg(17„/0.25)"-^ > 
0.62 at 95% c.L, and crg(ri^/0.25)°-5 > 0.60 at 99% c.L 
These constraints are for the largest sample, Mr < —22.3; 
we tried the analysis for the smaller subsamples in case 
splitting the sample to reach a higher mass threshold 
would give superior results, but unfortunately this was 
not the case. These constraints can be easily imple- 
mented in Markov Chain Monte Carlo codes when doing 
global parameter estimates in the form of a soft boundary 
on this parameter combination. For the median value of 
CTg for f2„i = 0.25, the minimum mass thresholds for the 
Mr < -22.3, < -22.6, and < -23 samples are 3.3, 6.9, 
and 14.1 x IO^^H^^Mq, respectively. 

The sources of error for the galaxy-galaxy weak lens- 
ing analysis have been discussed in the previous section 
and are described in more detail in [30[, and the overall 
calibration uncertainty of 8% has been applied to obtain 
the final results. Similarly, the errors associated with the 
cluster abundance are described in the previous section, 
and while we apply a correction for the clusters missed 
by the LRG selection, this is correction justified by our 
lensing analysis of Main sample luminous galaxies. One 
source of error that we likely do not need to worry about 
is the sampling variance, given that we work with a vol- 
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FIG. 4: The distribution of best-fit erg values in the bootstrap 
subsamples assuming fixed flm = 0.25, for the three luminos- 
ity threshold samples: Mr < —22.3 (black, solid); < —22.6 
(red, dotted); < —23 (blue, dashed). Each value should be 
interpreted as a lower limit due to the failure of modeling as- 
sumptions. We see that the brightest sample gives the weakest 
limits, suggesting that luminosity fails to trace the halo mass 
with low scatter at the bright end. 

ume of about a cubic Gpc with of order 10,000 clusters. 
For simplicity, we quote the constraint at z = 0, but the 
actual constraint is at z = 0.24, and for a given value of 
rim, we have translated between the two using a AC DM 
cosmological model. To translate to a different model, for 
example one with dark energy equation of state w differ- 
ent from —1, one must multiply by the ratio of the growth 
factors between these two redshifts for the two models. 
Over the range of w allowed by the current constraints 
and the range of redshifts of concern, these effects are 
well below the current level of statistical precision. 



IV. DISCUSSION 

We have presented a method that can robustly deter- 
mine a lower limit to the amplitude of fluctuations in the 
universe. The method consists of a statistical weak lens- 
ing analysis on a stacked cluster sample compared with 
theoretical predictions, assuming that the sample occu- 
pies the centers of all of the most massive halos for that 
model. The method does not attempt to establish a mass 
estimate for each cluster, and is thus significantly less ob- 
servationally demanding than the traditional methods of 
calibration based on a large cluster subsample for which 
the mass-observable relation is determined. As a result, 
only an inequality can be established. On the other hand. 
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unlike other methods, this method uses mass information 
extracted from lensing on the complete cluster sample, 
removing one of the major uncertainties in the traditional 
approach. Furthermore, the main difficulties in the tra- 
ditional analyses (knowledge of incompleteness and scat- 
ter in the mass-observable relation) are removed when 
attempting to place only a lower limit. We have applied 
the method to a sample of 16,000 cluster candidates from 
SDSS. We find as{fl,n/0.25)°-^ > 0.68 at 50% confidence 
level and as{n,j0.25f-^ > 0.62 at 95% c.L, where the 
error is entirely due to statistical and systematic errors 
in the lensing analysis, while all modeling uncertainties 
have been absorbed by the inequality. 

At first, it appears surprising that a cluster abundance 
analysis gives a robust lower limit. One effect of scatter 
in the mass-observable relation relative to no scatter is 
to increase the number of clusters at a given observable 
threshold, since for a symmetric error distribution around 
the mean, there are always more low mass clusters that 
scatter into the sample than high mass clusters that scat- 
ter out of the sample, a consequence of the steepness of 
the halo mass function. This increase in abundance is the 
so-called Eddington bias. Ignoring the scatter, one would 
therefore conclude that the derived value of cg is higher 
than the true value. Consequently, one sometimes finds 
statements that, ignoring scatter, the cluster abundance 
gives an upper limit on erg [111, EOl ■ 

However, such statements are only valid if the rela- 
tion between the mass and observable is extracted from 
a complete (or at least unbiased) sample of halos at a 
given halo mass threshold. In practice, for a given sam- 
ple (whether flux- or volume-limited), this goal is impos- 
sible because of scatter: we must define the sample used 
to determine the mass observable-relation with a given 
observable threshold. Consequently, low mass halos for 
which the observable is above the mean will be scattered 
in, and high mass halos for which the observable is below 
the mean will be scattered out. Since there are more of 
the former than of the latter due to the declining mass 
function, the derived mean halo mass at a given observ- 
able will be underestimated. This is the Malmquist bias 
effect that acts in the opposite direction to the Edding- 
ton bias, and leads to an underestimate of ag [4l'|. Note 
that Malmquist bias is present for both flux-limited and 
volume-limited samples, since the latter are always de- 
rived from the former using redshift information. While 
the two effects oppose each other, our results suggest 
that the Malmquist bias always dominates, so that only 
a robust lower limit can be established. We present an 
analytical derivation of this result in Appendix [Xj This 
result should be valid for other cluster abundance anal- 
yses where the effects of scatter are not explicitly taken 
into account, as long as the halo mass determination is 
reliable (and as long as the sample itself is used to de- 
termine the mass-observable relation, rather than sim- 
ulations). In practice, mass determination may not be 
reliable for those mass estimates that are not based on 
lensing, so one cannot conclude that all of the derived 



limits from cluster abundance with no attempt to correct 
for scatter should be interpreted as robust lower limits. 
Nonetheless, it is clear that scatter is a serious problem, 
and that it can lead to the opposite effect on ug from 
what is often assumed. 

Our lower limit on ag is relatively weak, even in the 
context of the WMAP3 results, which suggest a low am- 
plitude of fluctuations, cts ~ 0.75 [42], compared to previ- 
ous analyses giving 0.9 [43|. The cosmological constraints 
we find here are likely to be improved further with better 
cluster samples. First, the errors in the lensing analysis 
are significant, and the 95% limit is 0.06 lower in erg than 
the median 50% value, so a larger cluster sample or a 
deeper survey with more background galaxies would re- 
duce the statistical error. Second, in our analysis we 
work with relatively low mass clusters around IO^^Mq 
and densities between 10~^ and 10^^(/i/Mpc)"^. While 
the mass function is sensitive to the amplitude of fluc- 
tuations in this range, the sensitivity increases as one 
moves to higher mass clusters with lower number densi- 
ties, where the mass function exponentially decreases, as 
shown in figure [3] Third, using isolated LRG luminos- 
ity as a proxy for cluster mass is rather unsophisticated. 
While a relation between the halo mass and LRG lumi- 
nosity has been established over the range of masses be- 
tween 5 X 10^^ Mq and 2 x 10^"* Mq, it breaks down above 
that , so one cannot use this method to identify a sam- 
ple of very massive clusters in SDSS. While the scatter 
between LRG luminosity and halo mass has not been es- 
tablished, it is likely to be significant. Consequently, it 
is unlikely that the value of the amplitude of fluctuations 
derived in this paper is very close to the true value. In 
Appendix \X\ we derive the amplitude of the effect assum- 
ing a simple lognormal scatter. The bottom panel of flg- 
ureOshows that for a 50% scatter in the mass-observable 
relation, the underestimate of erg is 6-7%, relatively in- 
dependent of the mass above 10^'^h~^MQ. However, the 
scatter can be even larger and may not follow the simple 
lognormal model. For example, if we apply corrections 
to the observed lensing signal consistent with the results 
for the brightest luminosity sample in [2l| to account for 
the effects of scatter in M(L), the best-fit values of erg 
changes by 15%. We do not suggest that the model for 
M{L) scatter used in that work is a close enough ap- 
proximation to the real one that this correction should 
be applied; rather, we include this information to make it 
very clear that scatter is an important issue for the SDSS 
LRG sample that can seriously affect measured cosmo- 
logical parameters. 

In the future, it would be worth repeating this anal- 
ysis on samples of clusters with an even lower number 
density, for which the abundance is more sensitive to the 
amplitude of fluctuations. While the number of clusters 
will be lower, their lensing signal will be higher, so as 
long as the signal to noise remains high, such analysis 
will yield useful constraints that may improve upon the 
ones established in this paper. Many sophisticated algo- 
rithms have been developed to select clusters and deter- 
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mine their mass based on their richness of red galaxies 
with photometric redshifts [44]. A particularly promis- 
ing method in the context of the SDSS is MaxBCG [1], 
which identifies clusters by their concentration of galax- 
ies along the red galaxy color-redshift relation and estab- 
lishes the mass-richness relation across a broader range 
of halo mass. We expect that applying our method to 
this and other samples of clusters will yield lower lim- 
its on the amplitude that may be a useful complement 
to the traditional analyses of cluster abundances which 
yield very small formal errors, but may contain hidden 
systematic errors. 

Another possibility is to apply this method to the sam- 
ple of regular galaxies, tracing halo masses around 
10^^ Mq. These galaxies are well below the nonlinear 
mass, in the regime where the halo mass function is sen- 
sitive mostly to rimines + 3), where is the effective 
slope of the power spectrum at megaparsec scales (45| . 
Here again only a lower limit to this combination can 
be established. However, the errors on these parameters 
from other probes are already quite small, and it is not 
clear that such analysis would be competitive given the 
current errors in the lensing analysis. 

Is it possible that ultimately the lower limit we es- 



tablished here is all that will be achievable with cluster 
abundance studies, even when a reliable mass tracer such 
as weak lensing is used? This would be a major retreat 
compared to the expectations for future surveys [H, 0, 0] ■ 
In principle this option cannot be discarded, since there 
is always a possibility that a subset of dark matter halos 
is not detected with X-ray, optical or SZ surveys. More 
generally, the conditional probability distribution of halo 
mass at a given luminosity may be complicated and dif- 
ficult to establish using subsamples of data. While it is 
likely that further studies, both theoretical and observa- 
tional, will improve our knowledge, the ultimate limita- 
tions of the method are difficult to establish. It is there- 
fore important to know that at least some cosmological 
information extracted from cluster abundance is robust, 
even if it is only in a form of inequality. 
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APPENDIX A: ANALYTIC COMPARISON OF 
EDDINGTON VERSUS MALMQUIST BIAS 

Here, we analytically compare the magnitudes of Ed- 
dington and Malmquist biases with minimal assumptions 
for an analysis that proceeds as follows: determining 
the mass-observable relation from the lensing signal for 
the complete sample, and using the observed abimdancc 
with that mass normalization to constrain erg. Due to 
our strict cut in the observable, the steep decline of the 
mass function is going to (a) bias the average mass low 
(Malmquist bias) and (b) bias the abundance high (Ed- 
dington bias), whether the sample is flux- or volume- 
limited. We will now analytically prove our claim that 
the former effect dominates over the latter, leading to an 
underestimate of cts. We demonstrate this effect assum- 
ing a lognormal scatter in the mass-observable relation 
without placing a requirement on the size of the scatter. 

The parameters of our analysis are the true cluster 
mass M and the estimate of the mass M derived from 
a cluster observable such as richness, X-ray luminosity, 
etc. Due to onv use of a lognormal error distribution 
(i.e., constant relative rather than absolute scatter in the 
observable as a function of mass), our calculation uses the 
variables x = In (M/Mg) and x = IniM/Mo) for some 
arbitrary Mq. The relation between x and x is then 



p{x\x) 



1 



V2 



■.cxp[-{x-xf/2a% (Al) 



We assume the halo mass function can be expressed as 
a summation over power laws, so we perform the calcu- 
lations for a single power law and determine over what 
ranges of power law slope the calculation is valid. When 
using 



dn riQ / M 
we then find dn/dx = no exp [(a -|- l)x]. 



(A2) 



1. Observed abundance calculation 

First, we assume that the mass threshold is placed us- 
ing the observable by requiring x > Xc for some Xc- The 
observed number density is then 



n{x > Xc 



/•OO PC 



dx^p(x\x) = [ 



,-dn 
da;— 
da; 



(A3) 

where dh/dx is the observable mass function including 
effects of scatter. The integral over x can be performed 
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by completing the square to obtain 

— = no exp [(a + l)i] exp [cr^ (a + 1)^/2] 



dn{x) 
dx 



exp [cr2(a + 1)2/2] 



(A4) 



The latter term represents an increase in abundance due 
to scatter relative to the true halo mass function. Per- 
forming the integral over x in equation IA3I for a < — 1 
(true for all mass ranges of interest in this paper) yields 



n{x > Xc) 



exp [{a + l)x,] exp + 1)2/2] 

(A5) 



-ia + l) 
= n(x > Xc) exp [0-2(0, + 1)2/2] 



The last factor in this equation represents the Eddington 
bias (which is always greater than unity) and shows that 
the observed abundance is higher than it would have been 
at the same threshold in the absence of scatter. This 
effect is stronger for a steeper mass function because, for 
a fixed scatter, there are even more lower mass halos to 
scatter into the sample when the mass function is steeper. 
If the scatter is ignored then this effect causes one to 
overestimate erg from the observed abundance. 



2. True mean mass 

We now compute the true mean mass of the sample, 
which must also be affected by the lower mass halos scat- 
tering into the sample. This is the mass that will be 
measured from the sample assuming the mass estimator 
is correct. Since we have assumed that the scatter in 
mass is lognormal, i.e. gaussian in x, we compute the 
average of x rather than of the mass itself. 



1 



n{x > Xc) Jx 



dn 

dx I x—p{x\x)dx. (A6) 
dx 



The integral over x can again be performed by completing 
the square, and after carrying out the integral over x we 
obtain 



1 



1 



+ {a + l)cr2 (i) + {a + 1)0-2. (^7) 



The first term (x) is the estimated mean mass assum- 
ing our mass estimator. It equals the mean true mass 
in the absence of scatter. The final term represents the 
Malmquist bias; because a + 1 < 0, scatter causes the es- 
timated mean mass to be biased high relative to the true 
mean mass. Thus, if one neglected scatter, one would 
conclude from this sample that the mass estimator is bi- 
ased and would apply a correction from equation IA7I to 
rectify it, even though we have modeled the scatter be- 
tween the estimated and true value of x as a gaussian 



with zero mean. In the absence of an additional external 
calibration this is the only way to calibrate the estimator. 
While we have assumed that the mass estimator is given 
for every cluster in the sample, a random subsample of 
the complete sample would lead to the same result. 
a. Effects of ignoring scatter 

We now consider the competition between the com- 
puted Malmquist and Eddington biases in an analysis 
that ignores the effects of scatter. From lensing, we see 
a mean mass given by equation IA7I and determine an 
underestimated lower mass cutoff for our sample. 



Xc+ {a + l)o2 < Xc 



(A8) 



Next, we compute the expected abundance for that lower 
mass cutoff assuming no scatter in the mass-observable 
relation, i.e. 



. ^ ^ ^ dn 

n[x > Xc) = / dx — 
Jx^ dx 



no 



■ exp [{a + 1)3 



-ia+l) 

= n(a; > a;c)exp [(a + l)2o-2] (A9) 

Finally, we compare the true observed abundance in- 
cluding scatter, eauation lASi against the expected abun- 
dance in our analysis without scatter, equation IA9I We 
find 

n{x > Xc) _ exp [0-2 (a 4- 1)2/2] 
n{x > Xc) exp [o2(q! -f 1)2] 

= exp[-o2(a + 1)2/2] < 1. (AlO) 

Thus, for an arbitrary size lognormal scatter in the 
mass-observable relation, the Malmquist bias dominates 
over the Eddington bias, leading us to underestimate o-g- 
This conclusion is valid when the mass function can be 
expressed locally (near Xc) as a power law in mass with 
slope steeper than —1, which is true for all masses of 
interest in this work and, indeed, in most typical clus- 
ter abundance analysis. As a demonstration, figure O 
shows a (the power-law index of dn/dM) for — 0.25, 
as = 0.75, z = 0.25, and Ug ~ 1.0 for the mass ranges of 
interest, for which a < — 2. The bottom panel also shows 
the expected bias in o-g as a function of lower mass cutoff 
for fixed a = 0.1, 0.3, and 0.5 (no variation with mass) 
given the a from the top panel and the change in pre- 
dicted abundance with o-g from figure [S] This bias arises 
due to the dominance of Malmquist bias leading to an 
apparent underestimate of number density, which is in- 
terpreted in an analysis without scatter as a suppression 
of o-g. 
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FIG. 5: Top panel: power law index of dn/dM for the mass 
range of interest for this work. Bottom panel: resulting bias 
in o"8 when ignoring scatter in the mass-observable relation at 
the 10%, 30%, and 50% level. 



